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Abstract 

We explore gauge fields - strings duality by means of the loop equations and the zigzag symme- 
try. The results are striking and incomplete. Striking — because we find that the string ansatz 
proposed in [T|] satisfies gauge theory Schwinger- Dyson equations precisely at the critical di- 
mension D cr = 4. Incomplete — since we get these results only in the WKB approximation and 
only for a special class of contours. The ways to go beyond these limitations and in particular 
the OPE for operators defined on the loop are also discussed. 
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1 Introduction 



Gauge fields - strings duality is an old and fundamental subject. It underwent a rapid and 
fascinating development in the last three years. In this duality color-electric flux lines of gauge 
theory are described as certain relativistic strings. It has been shown in that the "natural 
habitat" for these strings (and thus for the flux lines) is a five-dimensional curved space with 
the metric 

ds 2 = d(p 2 + a 2 ((p)dx 2 . (1.1) 

The function a 2 {(f)) (which represents the running string tension) must be of a special type. 
In order to properly describe the zigzag-invariant Wilson loop, it must have a horizon, where 
o- 2 {4>h) = 0, and an infinity where a 2 (<f)i) = oo. The precise form of a 2 ((f>) is determined from 
the condition of conformal invariance on the world sheet. The gauge fields - strings duality was 
formulated in |1J as an isomorphism between the closed string vertex operators and the gauge 
invariant operators of gauge theory. 

Another development M was related to the M = 4 SYM theory and the 3-branes in the type 
IIB strings. It has been shown there that some of the SYM correlation functions (describing 
absorption of the external particles by the 3-brane) can be calculated by the use of the classical 



supergravity in which the 3-brane is replaced by the metric (\1 . 1| ) . It was further conjectured 
H that the correspondence extends beyond the SUGRA approximation and the function a 2 ((f)) 
in this case must correspond to the AdS space (a 2 ((f)) oc e a< ^). 

In the subsequent papers H it was shown how to implement the above-mentioned iso- 
morphism in the Af = 4 SYM theory. After that the correspondence in this case has been 
confirmed in an almost infinite number of papers. 

But all is not well. There is no real understanding (beyond the heuristic arguments of 0-0) 
why and, more importantly, when the correspondence works. Ideally, one would like to check 
that the Schwinger- Dyson equations of the Yang-Mills theory can be obtained from the string 
representation. This would clarify which string theories must be used for the various gauge 
theories. It is important to stress in this respect that the help of the D-branes in answering the 
above question is not always available. There are reasons to believe that in the most interesting 



non-supersymmetric cases the D-brane interpretation of the general cx-model metric ( |1 . 1|) is not 



possible. In these cases the Schwinger-Dyson equations (formulated as loop equations) is our 
only tool. 

In the vast literature on the subject there have been occasional attempts to explore the loop 
equations, but no conclusive results have been reached so far. 

The reason, we believe, lies in some common misconceptions concerning loop equations. It is 
generally thought that the loop operator is singular and can be applied only to the regularized 
and non-universal Wilson loop. This is not necessarily so. In || it has already been sketched 
how to implement the loop equation for the renormalized loops. 

In the present paper we shall apply the loop operator to the string functional integral. 
The results are incomplete but quite stunning. We will be dealing only with the very special 
contours — wavy lines (suggested in [jl[) and only in the WKB approximation. This is obviously 
incomplete. We will find in this case that, first, the loop operator is well defined in any 
dimension (implying the zigzag symmetry of the string representation), and second, at the 
critical dimension D CT = 4 the loop equation is satisfied. This we find quite stunning. 



2 The loop equation and the zigzag symmetry 

The loop equation has been introduced in and further explored in || and |§ (see also some 
related developments in ||, ]10[). We will give its derivation now, which can be used for the 
renormalized Wilson loop. It will contain some additional elements to the old works. The 
Wilson loop is given by 

W[C) = ^TrPexp j A^dx^ . (2.1) 

The averaging in this formula is performed with the Yang-Mills action 

1 
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S = ji- I TrC (dx) 



{F^v is the Yang-Mills field strength). The idea of the loop equation is to find an operation in 
the loop space which, being applied to the LHS of (|2.1|) , will give the Yang-Mills equations of 
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motion at the RHS. To implement this idea, consider the second variational derivative of V4^[C] 
5 2 W 



jTrp(F„ A (i( )> ))^(i( s '))^ x »' te »))i A ( s )i„( s ') (2.2) 



(where P is the ordering along the contour). 

The main idea of the loop equation is to separate the first term, and to introduce the 
"loop Laplacian" which acting on VK[C] gives the Yang-Mills equation of motion. This can be 
achieved in several different ways. The problem to overcome is the singularity of the second 
term at s — s', which makes it difficult to distinguish it from the first one. 

The "brute force" method would be to regularize the gauge theory with some cut-off A 
and the consider \s — s'\ <C 1/A. The second term is regular in this case, while the first one 
contains the 5-function which is easy to pick up. The price to pay is the non-universal arbitrary 
regularization. We cannot be sure that on the string side the regularization is the same and 
thus the comparison of gauge fields and strings becomes strictly speaking impossible. It has 
already been noticed in || that the way out of this difficult is to use the operator product 
expansion as s — > s'. 

Consider a set of operators {O n (x)} which are not color singlets. We can in general examine 
a gauge invariant amplitude 



G ni ...n N (Si, • • • S N ) = (TlP (O ni (x(Si))O n2 (x(s 2 )) • • • €~ 



Let us assume that the contour C is smooth and non-selfmtersecting. In this case we expect 
the OPE to have the formF] 

C m (x(s)) 

O n Msi))OnMs2)) =Y, Tr\ 7,XL 2 -A m °Ms)) (2.3) 

\x{Si) — X{S2)\ 1 2 

(where s = Sl ~^ S2 and A n are anomalous dimensions of the corresponding operators). Of course, 
in asymptotically free theories there are also powers of logarithms log |a?(si) — x{s2)\ in these 
formulas, which we do not display. 

In this way the OPE in the physical 4d space is transplanted to the one-dimensional contour 
C. The "open string" amplitudes G ni _ nN (sx, . . . sn) have power singularities at the coinciding 
1 Such operator products were considered before in [^6). 



points. There is no room for the contact terms proportional to 5(s — s') if we use renormalized 
correlators. Moreover, the possible contact terms in the x-space do not give contact terms in 
the s-space. Consider as an example a singularity 

5(x(s) - x(s')) = lim (4 - A) 1 . 

A^4 |x(s) — X(s')\ a 

As it was explained in ||, when one uses test functions, 

/ - *< s '» w ds ' = B< 4 - A ' / = ■ 

provided that the contour has no self-intersections, so that the only singularity in the integral 
comes from s = s'. It is also assumed that the integral is defined by analytic continuation. 

As a result of this discussion we conclude that the second variational derivative of the 
renormalized (and analytically regularized) Wilson loop has the form 

r /T^ WTO^O+E , ffiJ^W , (2.4) 

dx^sjdx^s') \x(s) — x(s')| 4 ^™ 

where C n (s,{x(s)}) oc <TrP (O n (x(s))ef Adx )). Since \x(s) - x(s')\ ss ^x 2 (s)|s - the 
nonlocal power-like behavior of the second term for s — > s' can be separated from the local 
s — s' singularity of the first one. The loop Laplacian is defined as the coefficient L(s)W in 
front of the 5-function, and the loop equation for non-selfintersecting contours has the form 

L{s)W = . (2.5) 



In our applications it will be convenient to write (|2.4j ) in the momentum space: 



This is an important formula. It shows that the Wilson loop belongs to a very special class 
of functionals, for which the asymptotics in ( [2.6|) does not contain integer even powers of p, 
corresponding to the derivatives of the 5-function. If it does, this means that the functional 
W[C7] cannot be represented as an ordered exponential ( |2.1| ). 

Some danger for this analysis would be presented by odd integers A n in ( |2.4|) , which would 
give a p k (logp + const) contribution. However, the lowest operators which can appear in ( |2.3| ), 
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even when protected by non-renormalization theorem, have even dimensions. An exception 
from this is the operator 'W a F^ which can appear in the OPE: 

F, x {x{s))F, a {x{s')) ~ ~* {S,) ]»Zf Xa ■ 

\x(s) — x(s')\ 

(We assume here the scenario in which the operator has the normal dimension as a result 
of some non-renormalization theorem.) However, the singular contribution from this term to 
( |2.2| ) vanishes after contraction with x\(s)x a (s'). 

The above analysis can be extended to the functionals represented by the iterated integrals 

W[C] = J2 [ T^(x(8 1 ) 1 ...x(8 n ))x^8 1 )...X^(8 n )d8 1 ...^ (2.7) 

n J S!<S2...<S n 

(of which ( |2.1D is a special case). Such functionals behave well under the application of the 
loop Laplacian (see |§). Hence, the second variational derivative for them has the structure 
(|2l 



Another feature of the functionals ( |2.7|) is the zigzag symmetry [![]. When we change s =>- 
a(s), they are invariant even if a'(s) changes sign (and is not a diffeomorphism) . It is likely 
that any zigzag-invariant functional can be represented in the form fl2.7|) , if it satisfies some 
extra analyticity requirements. Without them many other forms are possible, e.g. 

Ob qOC ^ Ob j^Ob 



L af3 



ds 



x 2 



At present the precise form of these requirements is not known, although a possibility is that 
the functional must be nonsingular with respect to x 2 . In the following, we will use the term 
"zigzag symmetry" referring to the form ( |2.7|) . According to the above discussion, there is a 
practical way for testing the zigzag symmetry in this sense: take the asymptotics ( j2.6| ) and 
check for the terms a p n with n — positive even integer. If they are present, there is no place 
for the zigzag symmetry, and the functional does not represent any Yang-Mills theory. This 
conclusion, together with the formula ( |2.6| ), will be extensively used below. 



3 The string representation: D-branes and/or a- model. 

Let us summarize the facts about string representation of the Wilson loop. According to [|TJ 
we must consider a 5d background with the metric ( |1.1|) . The origin of the fifth dimension is 
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the Liouville degree of freedom which is unavoidable in the non-critical strings. The cr-model 
action of this string is given by 

S = ljd 2 ^g a \0GMN(z(0)d a z M d b z N + ^z(0) (2) R(9)V9 

+ e ab B M N{z{C))daZ M dbZ N + other background fields + fermions. (3.1) 

We introduced the 5d variable z M = We are also considering a fermionic non-critical 

string. The metric Gmn has the form obtained by the obvious change of variables in ( |1.1|) : 

ds 2 = p(y)(dy 2 + dx 2 ) . (3.2) 



As usual, the background fields are determined from the /3-function equation [|11J| , expressing 
the g ab -independence of the theory. 

It is important to have the NSR fermions for the gauge fields - strings duality. It is also 
necessary to perform the non-chiral GSO-projection in the functional integral. The whole 



matter has been treated in the paper |12|]. Here we will just remind the reader the underlying 
logic with a few additional comments. 

There are two types of objects we would like to consider. First, there are the closed strings 
amplitudes. We introduce the closed string vertex operators, e.g. 

= [ d\ U^y{£))d a x»d a A e**® 



and claim that for the proper choice of the metric p(y) the correlators of V's on the string side 
are equal to the correlators of the gauge invariant operators on the Yang-Mills side. In other 
words, we conjecture an isomorphism 

Closed string Gauge invariant 

<*=5> (3.3) 
states operators. 

The intuitive reason for this conjecture is that gauge invariant operators like TtF 2 u can be 
associated with tiny closed flux lines, represented by small Wilson loops. Indeed, we create a 
flux line by acting on the vacuum with the operator 



TrP exp ^Adx^ Q I + C^ Xa Tr(F^F Xa ) + ... 

c 
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The propagation of these loops forms the closed world surfaces with punctures represented 



by the vertex operators. This means that if we find the right action ( |3.1| ), the closed string 
^-matrix will be equal to the correlation functions of the local gauge invariant operators. 

But how to discover the right action? This question does not have a complete answer 
at present. There are two approaches the problem. One approach applicable to the 

supersymmetric Yang-Mills theories, is to start with the 3-branes of the Type IIB superstring. 
Then, on the one hand, it is known that the low energy excitation of the stack of N three-branes 
are described by the supersymmetric Yang-Mills theory with the SU(N) gauge group. On the 
other hand, 3-branes can be presumably replaced by the supergravity background the generate, 
and thus one expects and checks the isomorphism ( |3.3| ) . 

The main problem with this approach is that it is far from being clear that non-supersymmetric 
theories can be described in the D-brane language. Also, even in the supersymmetric case, so far 



there is no real derivation of the isomorphism (|3.3|) from the first principles, although numerous 
checks have been performed. 

Another approach [|lj is to try to formulate the conditions under which the cx-model satisfies 
the Schwinger-Dyson equations of the Yang-Mills theory. The basic object in this approach is 
the Wilson loop, defined by 

W[x(s)\ = j Vy(OVx(Oe- s[m ' ym . 

AdD = x(s) 
V\dD =? 

Here we take the world sheet with a disc topology (appropriate for large N). At the boundary 
of the disc we fixed x by the condition x\qd = x(s). Now, we have to determine y\dD and the 
background fields so that the loop equations are satisfied. The basic idea of was to notice 
that the metric p(y) must have horizon and/or infinity points, so that 

p(y H ) = and p(yj) = oo . 

These are the only positions where we can place the contour. Indeed, if we choose some other 
y\dD, we will encounter the lack of the zigzag symmetry coming from the fact that the first 
term in (]7]) is (due to the presence of y/g) invariant only under diffeomorphisms, but not under 
the zigzag transformations. Hence it must be either zero or infinite at the boundary. A more 



precise condition was given in [O. It was stated there that in order to describe gauge theory 
we must find the metric p(y), the other background fields, and the boundary conditions y\dD 
such that the vector vertex operators, defined at the boundary, form a closed algebra. Physically 
this means that open strings in this theory correspond to vector gluons. Also, the vector vertex 
operators are selected by the zigzag symmetry, since they are the only ones defined without the 
world sheet metric. 

When the D-brane description is applicable, the above condition is satisfied. But it is 
important that the condition be formulated without any reference to D-branes. 

In the case of the M = 4 SYM the D-branes and conformal symmetry considerations lead 
to the formula [|] 

p(y) = \[qJ^n ■ 1 . (3.4) 

It has been shown in [f|], [||] for the local operators and in f[3 |, |fL4 | for the loops that the natural 
location for the gauge theory quantities is the infinity in this AdS space, that is y — > 0. 

In this paper we shall adopt this prescription, although there is no real proof for it. In a 
sense, our check of the loop equations in this paper can be considered as such (incomplete) 
proof. 

One should keep in mind, however, that there is another zigzag-symmetric location, y — oo. 
This point has no geometrical significance in the Euclidean signature (the Lobachevsky space) 
but is meaningful in the Lorentz signature (that is in the AdS space). In the latter case this 
point represents a horizon behind which the D-branes are hidden. This horizon is physically 
important because the absorption by the D-branes can be accounted for as the disappearance 
of the closed strings behind the horizon 0. It is therefore natural that the Wilson loop C on 
the D-brane can be represented as a Wilson loop at the horizon []]J. In this case, however, the 
loop generates the .B^-field with B^ v = B flv (x, {C}) |T[. At present we do not know how to 
determine this field. We hope that the two types of boundary conditions yield the same result 
and will be studying the case y — 0. 

The string representation involves various background fields. For example, since we need 
NSR fermions to exclude the boundary tachion (violating the zigzag symmetry), we have RR 
fields in our closed string spectrum, and they form a condensate in order to stabilize the AdS 
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space. The concrete types and forms of background fields depend on the type of the gauge 
theory we are working with. The problem of finding the precise background in each particular 
case is not yet completely solved. Of course to check the complete loop equations we have to 
solve the above problem first. However, in the conformal cases it is possible to retreat to the 
quasiclassical domain by taking gy M N 3> 1. In this case the string action takes the form: 

Smivm = J {(d a x) 2 + (d a y) 2 ) +0(1) , (3.5) 

where the terms 0(1) contain NSR fermions and RR fields. We arrive at the conclusion that 
in the WKB limit, the Wilson loop is given by JT3|, [|TJ 



W[C] OC e -V#M^min[C] _ (3 

Here 

A min [C] = min i J ^ {(d a x) 2 + (d a y) 2 ) (3.7) 

is the minimal area of a surface bounded by the loop C. In the following sections we will 
analyze the action of the loop operator on this functional. 



4 Minimal area in the Lobachevsky space 

Our contour is located at infinity (the absolute) of the Lobachevsky space. There are some 
general features of such minimal areas which we discuss in this section. The equations of motion 
for the action ( |3.7|) have the form: 

d a [\d a x] = , 

W J (4.1) 

d 2 y = - {(d a y) 2 - (d a x) 2 ) . 
y 

We also have to impose the Virasoro constraints 

(<9i£) 2 + (d lV ) 2 = (d 2 x) 2 + (d 2 y) 2 , ^ 
d\x d 2 x + d\y d<iy = . 

We are looking for the solutions satisfying the conditions 

x(a, 0) = c(ct) , 
|/M) = 0. 
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(We renamed the variables: £ x = a, £ 2 = r.) It is easy to see that the expansion in r has the 
form: 



* = 3» + \f{e)r 2 + \g(°V + ■■■ 
y = a {a)r + h{a)r 3 + . . . 

After substituting this expansion into ( |4. 1| ) and (O) we obtain 



(4.3) 



a 2 (a) 



/ dc \ 2 d ( d a c 



m Hz z 



Vrfay da V(<9 ff c) 2 
This guarantees that the leading term in the energy-momentum tensor 



(4.4) 



vanishes; we have 



0i 



(d T xd a x + d T yd a y) 



a?T 



fc + aa! + — (go) + . . . 
J cr 



and due to (fOj) the first term vanishes, while the second one gives 

1 



0i 



(*0 • 



Analogously: 



flu. = [(^ f ) 2 + ( d ry) 2 - ( d ^) 2 - ( 9 «y) 2 ] 



/ 2 + 2a6 - c'f ' - a' 2 



(4.5) 



The action ( p.7|) calculated on the classical solution (|4.3|) has the following structure (encoun- 



tered previously in [ 13 1 in a special case) 



A min [C] = ^ + A[c(a)] , 



(4.6) 



where we introduced the cut-off y min = a(a)r min = e. In this formula L[C] is the length of the 
contour C and ^4[c(a)] is a finite functional. 

Our main interest is to derive variational equations for ^4[C]. One might think that the 
usual Hamilton- Jacobi equations following from the conditions 0jji = 8±± = will give us a 
closed equation for A. Unfortunately, life is not so simple. The above would be the case is we 
were able to find the coefficients g and b in terms of c. Indeed, it is easy to see that 

5A g{a) 



5c(a) 



a? (a) 
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and if b were known, the equation for A would follow from (|4.5|) . 

There is an unpleasant surprise, however. Further iterations of (|4.1|) reveal that the functions 
g(a) and b(a) in ( f4.3|) are not determined by the small r-expansion and can be kept arbitrary. 
They are fixed by the global condition for the absence of singularities at finite r, and thus it 
is hard to find them explicitly. Because of this difficulty we will choose an alternative way of 
finding ^4[C]. Namely, we will consider a special type of contours — wavy lines, and will develop 
a method of successive approximations for ^4[C]. We will also see another derivation of (|4.6| ). 



5 The theory of wavy lines 

It was already suggested in |L] that it is instructive to consider wavy lines instead of general 
contours, namely to look at a curve 

x 1 (s)=s, x*(s) = 0*(s), i = 2,...D, 

— * 

and to assume that <fi l {s) are small. Below we shall find the expansion of ^4.[0(s)] up to the 
fourth order and then explore the action of the loop Laplacian. To perform this expansion it is 
convenient to consider the standard Hamilton- Jacobi equation for the minimal surface which 
has the well-known general form: 

rMN( ,JA_ 8A dz M dz N 
G [Z) 6z M (s)6z"(s) ~^ MN ^blb- 

For the Poincare metric this equation becomes 



> fcwy ,5.1) 



$y( s ) J \Sx(s) J y A (s) I \ds J \ds 



SA 



We want to explore the limit y(s) = y — > 0. We can do this by solving (|Q|) with respect to s 
and by noticing that 



OA f , 6A 
ds- 



dy J 5y(s) 



(5.2) 

y(s)=y 



We have from ( 5.1|) and (|5.2| 



dA If, I fdx\ 2 . ( SA 
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We see directly from ( [5.3|) that A(y) behaves like 



A(y) » JM + (1). 
To get further information we have to look at the following expansion of A: 
A = S ^-] I ds 1 ...ds n T h _ in (si,...s n \y){<p ll (s 1 )...<p in (s n )) 

n 

Expansion of (]5.3|) together with the reparametrization invariance equation 



dx^ 5 A 



5A 



ds Sx^s) 5xi(s 
will give us equations for the T's. Let us expand Q5.3Q : 



dcj) 6 A 
xi=s as Sd) 



(5.4) 



OA 



Lc, 



dy y 2 + 2 J y 71 y 2 ^ ) dS 



1 -,2 



(5.5) 



where n = 5 A/ 5(f). If we perform the Fourier transform in s, we obtain the following equations 
for the coefficient functions: 

,2 



f dT 2 _ 2 2 _ r 

dy y 2 

' 4 



dv 4 



J^fe) r 4(Pl,---P4) ~ B 4 (p 1} . . . p 4 ) . 



(5.6) 



We assume here that 



A =7 + ^ / r 2 (p)(^- P ) ^ 
1 



(5.7) 



r 4 (pi,...P4) U> 

) 5 (X^*) dpi . . .dp A + . . . 
The only slightly complicated structure is the polynomial B 4 which is found from (|5.5| ): 

(2n)B 4 (p u ...p 4 ) = C 4 (pi, ...p 4 ) + D(p u p 2 ) + D(p 3 ,p 4 ) - 

-D(p u p 3 ) - D(p u p 4 ) - D(p 2 ,p 3 ) - D(p 2 ,p 3 ) . (5.8) 



In this formula 



C4 = {P\P2P?,Pa) —Z + V 6 (Wi^ 2 ^3^4) 

yi 

D(pi,p 2 ) = (PiP2u 3 u 4 )y 2 , 



1.3 



where we introduced the notation u>i = T 2 (pi). 

The equations ( |5.6|) are easy to solve. From the first one we get 



p 2 

r 2 (p, y) = u(p, y) = 5.9 
2/(1 + my) 

(this solution is the only one which is positive and regular for y — > oo.) From the second one, 

POO 

T 4 = / dye-f° d ™ 2 ^)B 4 ( Pl ,...p 4 \y) = 
Jy 

POO 

= / dy 11(1+^)6-^^13^..^^). (5.10) 

Jy i 

We see that T 4 also has the structure of ( |5.8|) : 

(27T)r 4 = F(pi, ...Pi) + $12 + $34 - $13 - $14 - $23 ~ $24- (5.11) 



Taking the integral (5.10) and separating the finite part for y — > gives 



2 2 2 2 
PlP 2 »> 



p _ f 2 ^2^4 + ^ | ei^Q Ay^ 1 \ | J2i<j\Pi\-\Pj\ A \ 

V A3 ^ 2 V \Pi\J ^P\PlPzP<i PlP2P2,Pi) 

(A = ^ bill e i — s § n This completes the calculation of r 4 (pi, . . .p 4 ). 

Let us notice that an alternative way to obtain these formulas is to use the Monge gauge in 
the expression for the minimal area: X\ = a, y = r, 0j = 0j(cr, r), 



A = J tlLyJ ! + $* + $* + 0202 _ j^) 2 . (5.13) 
The linear equation for <p, 



o 1 r I ' o cri 



has the solution 

&(P,T) = (\p\rf 2 K 3/2 (\p\T)$(p) = (1 + |p|r)e-l^0(p). 



Expanding (|5.13 ), we arrive once again at the expression ( 5. ICQ . However, the Hamilton- Jacobi 



method has some general advantages. 
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6 The second derivative of the minimal area 



After the divergent part of (|4.6|) is absorbed into the mass renormalization of the test particle, 
we are left with the functional 



W[C\ = e-V&^" A[Ci . (6.1) 



Its second variational derivative has the form 

5 2 W ( , „ SA 5A nr—, S 2 A 



2 Ar dA dA rz — ~ 



w . 



The first term in brackets has no singularity for s — > s'. This fact was considered in as 
a check that the loop equation is satisfied. Notice however that in this order an arbitrary 
functional A will pass this check. That is why it is necessary to consider the second term, 
which we proceed to calculate. 

The full second variational derivative of A consists of the longitudinal and transverse parts: 

6 2 A = 6 2 A 5 2 A 
Sx^Sx^s') Sxi(s)Sx 1 (s') 5<p{s)5<p{s') ' 



The transverse part (in momentum representation) can be read off directly from (|5.7|) . For the 
term quadratic in (p we have 

8 2 A^ 

lim -j- —j- r = (D- l)8(q)T 2 (p) = (1 - D)5(q)\p\ 3 , 

where T 2 (p) = — \p\ 3 is the finite part of T 2 (p, y). 

Let us now consider the quartic term. The special structure in ([5.11 ) leads to the formula 



(2?r) = -- / H(k,k',pi,p 2 ) \(f) Pl (f) P2 ) S(k + k' + pi + p 2 ) dp x dp 2 , 

d(p{k)0(j){k ) A J v 7 

H = (D + l)F(k, k',p 1: p 2 ) + (D - 3) k') + ${p u p 2 )} 

-(D - 1) [$(fe,pi) + ®(k,p 2 ) + $(k',pi) + $(k',p 2 )] . (6.3) 



Let us first analyze the case k = —k', which corresponds to taking q = in fl2.6|) . In this case 
the third term in the expression for H in (|6.3j ) drops out, because $(/c,p) is antisymmetric 
with respect to k — > —k. The remaining terms have to be expanded for k — > oo. Because of 



15 



homogeneity, it makes sense to express them as functions of x = k/p, where p\ = —p2 = p, 
which amounts to taking p = 1 in fl5.12| ). We have: 

F(k,-k, P ,- P ) = M 5 {^ + ^ + x2) - + «i 

= H ={-^-x + (I)}, 
= bl 5 {-^ 3 -x + 2 + oQ } . 

We substitute this into (|6.3|) and obtain the desired asymptotic expansion in k — » oo: 



8<f){k)8<f>{-k) 



2tt 



+ (3-D)|p| 5 }(0 p 0_ p ) * + • (M 

In general, when we do not assume that q = in ( |2.6| ), we have to take all the terms in (|6.3| ) 
into account. The corresponding formula for the asymptotic expansion (of which ( |6.4| ) is a 
partial case) can be obtained by means of straightforward but lengthy calculations and has the 
form: 

5(f) (f + k) 8$ (| - k) 2tt 7 I PlM 1 

+ (^-j^PipI + ^^-Plp^j \k\ + ((4 - D)pl\p 2 \ 3 + P!P 2 \P2?) } 

x (4> pi <j)p^ 5(pi + p 2 + g) rfpirfp2 + O • (6.5) 



To complete the calculation of ( |6.2| ), we have to compute the longitudinal part (the first 
term in ( |6.2j )). This part is not immediately visible from ( p.7| ) and has to be recovered by the 
use of the reparametrization invariance. The required identity can be easily derived from (|5.4| ) 
and has the form: 

' s )<m s Jtttttttt: -o{s- s )0(s)— 



8xi(s)8xi(s') ' 8(j)i(s)5(J)k(s') ds 

In the momentum representation: 

( 27r )"F — 77T7 777V = ~ / dp dp' pp'<f)i(p)<f) k (p ,x 



5xi (k) 5xi (A/) J * 0(pi(p + k)8(j)k(p' + k') 

f , ^ 8A 
I dpp{p + k + fc')0(p)— — . (6.6) 

J 5(b(v + k + k') 



6<j){p + k + k') 
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In the approximation we are working with, only A^ contributes to the RHS of (O). The 
result is: 

5 2 A 



1 

2^ 



Pi ~ P2 



k 



- \P2 



5x x (f + k) 5xi (f - k) 

XP1P2 (<Sn0 P2 ) 5(Pi +P2 + q) dp x dp 2 

=k^oo ^ J{\k\ 3 + \{Pi ~P2?\k\ -PiP 2 \P2\ 3 )x ■■■ 
The final expression for the asymptotics of the second variational derivative is thus: 



+ 



£) - 1 



3L>-9 



|fc| + (4- J D)pf|p 2 | 



x ( pi p2 ) 5(pi +p 2 + q) d Pl dp 2 + O ( ^ ] + 0( 1 5 



In the next section we will analyze this result. 



7 Interpretation and discussion of the result 



(6.7) 



(6.{ 



Our main result is contained in the formula ( |6.^ ). It shows, first of all, that the second 
variational derivative has the expected form fl2.4T) . Namely, dangerous terms oc 5"(s — s') 
(which would manifest themselves in (|6.8|) as terms cx k 2 ) cancel for all D. The presence of 
those terms would imply that our functional is not zigzag- invariant, that is not presentable in 



the form (277 ). 

The next result following from (|6.8| ) concerns the loop operator L q . Using (|2.6|) and picking 
up k° terms in 



we get 

4- D 



LgA 



2tt 



Pl\P2 



Pl (f) P2 J 5{pi + p 2 + q) dpxdp 2 



(7.1) 



This shows that at D = 4 the loop equation is satisfied (at least in our approximation)! We 
will discuss the significance of this fact in the next section. 

Now let us perform another test. Consider once again the OPE ( p.3|) and ( |2.4|) and let us try 



to determine its contribution to (6.S). If we assume that in conformal theory in the WKB limit 
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the field strength keeps its normal dimension 2 (this will be actually more of the conclusion 
than of the assumption), we get 

jp / \tp / \ n ^ Act i n ( Xl ~~ x 2)a(^i - x 2 ) a 

\X\ — X2\ \Xi — X 2 \° 

. n Oi - x 2 )(x(xx - x 2 ) ll F lur) (x) (xi - x 2 )^J^F Xa (x) 

+ ° 3 \z — rii + ° 4 1 12 • u-^> 

Fl — x 2| \X\ — x 2 \ 

(x = Xl 1 2 X2 ). The last two terms give no singular contribution to ( |2.2j ). The first two give after 
some calculations: 

5 2 W 1 d + C 2 + 1 / (Ci + C 2 )(xx) + + C 2 (xx) 2 



Sx^Sx^a') | s - s'| 4 x 2 \s - s'\ 2 V 12x 4 4i 4 4x 4 

The derivatives in the RHS are taken at the point s = In the wavy line approximation we 

get 

(7.3) 

We have to compare this behavior with our formula ( |6.8D . Picking up the terms oc |fc| 3 and \k\ 
and using the Fourier transform identities 

|s — s'| 4 6|s — 

to go back to the s-representation, we have 

^ K i (1 _ m _ ^ + i flz^ f + ^ . (, 4) 



Sx^Sx^s') |s-s'| 4V /v ^' \s-s'\ 2 \ 12 
We see that this formula can be put in complete agreement with ( |7.3|) by taking C\ = D — 3, 



C 2 = 2. The second term in ( [7.4| j may be modified if the theory contains a scalar operator of 
dimension 2. That can change the constants C\ and C 2 . At the same time, the structure of 
the first term cannot be modified by anything and provides a strong check for the consistency 
of our approach. Notice that it also predicts that the dimension of F^ is not renormalized. 



8 Conclusions and outlook 



The main efforts of this work were directed towards the development of new techniques for 
checking the loop equations in string theory We managed to apply the loop Laplacian to the 
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minimal area in the Lobachevsky space and to show for the first time that the equations of 
motion of gauge theory are satisfied by string theory, at least in the WKB approximation. 
This point perhaps requires some clarifications. Namely, we looked at the Wilson loop ( ]2.1p in 
conformal versions of the Yang-Mills theory. These versions unavoidably contain other fields. 
Hence we expect that 

V M F M „ = J v ± 

(where J u is the current generated by those fields). So, what is the meaning of finding that 
LA min = at D = 4? 

There are several possible interpretations of this fact. Let us begin with the unpleasant one 
(in which we do not believe). It may be that our result is just a fluke, while if we proceed to 
higher orders in "waviness" of our contour, the loop equations will not be satisfied. Further 
progress will be difficult in this case. 

Let us take another, optimistic view. Various gauge theories presumably have a string 
representation with the background (1.1). If the theory is conformal, ( |1 . 1| ) must describe the 



AdS space. What distinguishes various theories is not the metric but other background fields 
and also the field content on the world sheet. However, in the WKB limit gyM^ ~~ * 00 the 
asymptotics of the Wilson loop is given by the formula ( |3.6|) and is the same for all conformal 
theories. If this is the case, it means that the current in the above formula is negligible in the 
WKB limit, and our result actually checks the universal law ( |3.6| ). Notice also that we are 



considering the standard Wilson loop and not its modified version suggested in fT4]| . Once 
again, with the above philosophy this modification is irrelevant in the WKB limit. 

To verify the above assertions it is necessary to go beyond our wavy line approximation. We 
believe that this can be done by a more general treatment of the Hamilton- Jacobi equations. 
Conformal invariance of the loop Laplacian (which we discuss in the Appendix) should play 
an important role in this analysis. Alternatively, one can study the second variation of the 
functional fl3.5| ) by developing the short distance expansion of the Green functions for equations 
( |1 . 1| ) . It is conceivable that by this method it will be possible to relate OPE on the world sheet 
and OPE in gauge theory. 

This brings us to a more difficult problem of going beyond WKB approximation. Again, the 
method of wavy lines may be useful here, but at the moment we do not know how to evaluate 
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the action of L q on quantum corrections to our formula. 



In the case of non-conformal theories the metric has the form 



ds 2 = /(logy) ( 




The problem for our analysis is not so much the function / (it is easy to generalize our consider- 
ations to this case) as the absence of the WKB domain. It may be helpful to notice that instead 
of considering non-conformal case in 4d, one can, in the asymptotically free theories, shift to 
D = 4 + e, when these theories become conformal. Perhaps, the classical part of L q W oc e 
will be canceled by the quantum fluctuations and will provide us with the "Holy Grail" of this 
subject — the space-time /3-function. 
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Appendix: Conformal invariance of the loop equation 

The basic property of the functional A (and hence of W) is its conformal invariance: 



(where / is a conformal transformation, say f^{x) = x^/x 2 ). This is true for any D. Although 
this invariance is to be expected, it is not entirely obvious, since isometries of the Lobachevsky 
space (corresponding to conformal transformations on the boundary) will change the cut-off e in 
(|4.6| ). To check the invariance, we extend / to the (D + l)-dimensional Lobachevsky isometry: 



Consider the minimal surface M bounded by C. Since F is an isometry, F(M) is the corre- 
sponding surface for f{C). Moreover, 



A[C] = A[f(C)} 



(8.1) 




(8.2) 



Area[M e ] = Area[F(M e )] 



(8.3) 
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(where M e is the surface M cut off at the hight e). For small e, the LHS is equal to L[C]/e+*4.[C]. 
According to ( |3.2| ), the constant cut-off e is transformed by F to the variable cut-off e/x 2 . To 
calculate the RHS of flS.3| ), it is convenient to introduce an auxiliary constant cut-off e' on 
F(M), so that Area[F(M e )] becomes equal to Area.[F(M) e >] plus the area of a narrow strip 
between the variable and constant cut-offs. Now (|8.3|) implies: 



no], = wm +Aim]+ i *&w) is i 

e,e'^o e' Jc ds J e 



dy 



e/x(s) 2 V 2 



Calculating the integral, we see that the singular terms cancel, and we get ( |8.f| ). 

It is now natural to ask if the loop Laplacian transforms is conformally invariant, i.e. com- 
mutes with conformal transformations. This property can be written as the equality 

L{s)U f [C}=p(L{s)u)[f{C)} (8.4) 

valid for any (reparametrization invariant) functional U[C] (where Uf[C] = U[f(C)] is the 
functional U transformed by a conformal transformation; p is some factor). It turns out that 
(|8.4j ) is true if and only if D = 4. To prove this, consider the relation: 

62Uf[x{s)] - WW(^) »nxL + 52 /a(x(,)) 7 |^-5(, - s') . 



For conformal transformations we have 

Now, we must collect terms proportional to S(s — s'). This gives 
L(s)U f [C]5(s - s') = p(f) (L(s)u) [f(C)]5(s - s'] 



+ d 2 f x ^r6(s - s') + d,f x (x(sMJ a (x(s'))N [X(7] 5'( 



s — s 



where iV^] is the coefficient in front of the (^'-function in the second variational derivative: 



5f x (s)6f a (s>) 



From the condition of the reparametrization invariance s ^ s ^ = we have the identity 

SU 



5x\(s) 
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This gives the transformation law 

L(s)U f [C] = p(f) (l(s)u) [f(C)] + CVU/) N[>»] . 

= (d 2 f [x ) (d,f a] ) - (d a dj [x ) (d a f a] ) . (8.5) 

Substituting = x^/x 2 , we find that fi(/) oc (D — 4), and thus ( |8.4| ) is true if and only if 
D = 4. 

As a consequence of the above discussion, for D = 4 the loop equation for the Wilson loop 
(|6.1| ) is conformally invariant: 

L{s)W[f{C)\ = pL(s)W[C] . (8.6) 

It follows that in order to check the equation at a point x(s) of a contour C, we are allowed to 
first apply a conformal transformation, say with the purpose of simplifying the behavior of the 
contour at the point we are looking at. Although this observation does not play any significant 
role when working with wavy lines, it might become important for general contours. 

For D ^4 the presence of nonzero additional term in the RHS of the transformation law 
( |8.5| ) implies that the loop equation for ( |0| ) cannot be satisfied in this case. 



References 



[i] 

[2] 
[3] 
[4] 
[5] 
[6] 
[7] 



A. Polyakov, Proceedings of "Strings 97" , |hep-th/97lIUU2" 



I. Klebanov, Nucl. Phys. B496 (1997) 231, |hep-th/9702076 



J. Maldacena, Adv. Theor. Math. Phys. 2 (1998) 231, |hep-th/9711200 



S. Gubser, I. Klebanov, A. Polyakov, Phys. Lett. B428 (1998) 105, |hep-th/9802TU9" 
E. Witten, Adv. Theor. Math. Phys. 2 (1998) 253, |hep-th/ 9802150 
A. Polyakov, Nucl. Phys. B164 (1980) 171 
A. Polyakov, Phys. Lett. B82 (1979) 247 
Yu. Makeenko, A. Migdal, Phys. Lett. B88 (1979) 135 

22 



[9] 



Y. Nambu, Phys. Lett. B80 (1979) 372 



[10] 

[11] 
[12] 

[13] 

[14] 

[15] 



J. Gervais, A. Neveu, Phys. Lett. B80 (1979) 255 

C. Callan, E. Martinec, M. Perry, D. Friedlan, Nucl. Phys. B262 (1985) 593 



A. Polyakov, J. Mod. Phys. A14 (1999) 645, |hep-th/9809057 



J. Maldacena, Phys. Rev. Lett. 80 (1998) 4859, |hep-th/9803002 
S.-J. Rey, J. Yee, |hep-th/980300T 



N. Drukker, D.J. Gross, H. Ooguri, Phys. Rev. D60 (1999) 125006, |hep-th/9904H)T 



[16] 



H. Dorn, Fortschr. Phys. 34 (1986) 11 



23 



